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We consider the most general form ol extended Hubbard Hamiltonian conserving the total spin 
and number of electrons, and find all the 1-dimensional completely integrable models which can be 
derived from first degree polynomial solution of the Yang-Baxter equation. It is shown that such 
models are 96. They are identified with the 16-dimensional representation of a new class of solutions 
of symmetric group relations, acting as generalized permutators. As particular examples, the EKS 
and some other known models are obtained. 
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Recently the interest in low-dimensional correlated 
] electron systems has motivated some attention on the 
search for completely integrable 1-dimensional extended 
Hubbard models Q Q . The physical relevance of the 
latter has been realized in very different context, from 
high-Tj, materials to quantum wires Determining 
\ the integrability of model Hamiltonians can yield a wide 
set of rigorous results, very rare for interacting many- 
body systems, as the exact ground state phase diagram 
and the evaluation of correlation functions. In general, 
the complete integrability of a model is proved by pro- 
viding the underlying i?-matrix, solution of the Yang- 
Baxter equation (YBE), which is related to the factoriz- 
ability of the scattering matrix; whereas its spectrum can 
be obtained by diagonalizing the corresponding Trans- 
fer matrix. To achieve both the above purposes, the 
Quantum Inverse Scattering Method (QISM, see Q and 
references therein) has turned out to be a very power- 
ful tool. Within this scheme, the standard approach is 
to impose some appropriate (quantum) symmetries onto 
I the i?— matrix, and then to derive from it the integrable 
Hamiltonian sharing the same symmetries. Whether the 
model so obtained does actually correspond to a physical 
system or not, it is always an a posteriori check to do. 
On the contrary, when having a specific physically inter- 
esting Hamiltonian H as input, one may wonder whether 
a i?-matrix can be obtained as an output. To this pur- 
pose, in a recent paper |^ it was noticed how one can ex- 
plicitly build up the i?— matrix as a polynomial of degree 
p in the spectral parameter, with the requirement that its 
first derivative is precisely H . Constraints on the coeffi- 
cients of the power expansion were obtained, by imposing 
the YBE, in the standard form depending on the differ- 
ence of spectral parameters. In particular, it was shown 



that, in all cases in which the expansion was finite, the 
highest 3p-order equation is the YBE with no spectral 
parameter, i.e. the braid-group relation ||^, whereas the 
2p-order specializes this relation to the symmetric group 
case. 

Starting from this observation, in the present letter 
we consider the most general form of the (isotropic) ex- 
tended Hubbard model preserving the spin components 
and the number of electrons. By imposing that the lo- 
cal Hamiltonian satisfies the symmetric group relations 
we derive the constraints among the physical parameters 
which guarantee its integrability, the i?-matrix simply be- 
ing in this case a first-degree polynomial. In so doing we 
both obtain some known results, for instance the EKS 
model Hamiltonian, and prove the integrability of many 
other cases. It is realized that the solutions act as gener- 
alized permutators (see below) on the chain sites. At the 
same time, we find all the solutions of symmetric group 
relations having the form of a generalized permutator, 
and show that the above integrable extended Hubbard 
models are precisely those symmetric group generalized 
permutators (in the 16-dimensional representation) that 
preserve the Hamiltonian symmetries. 
Interestingly our results are obtained without imposing 
any extra symmetry constraints. Nevertheless, like in 
the quantum symmetric cases, they have a quite sim- 
ple structure, deriving from the underlying structure of 
generalized permutator solutions of the symmetric group 
relations. 

The most general form of the extended Hubbard model 
JlOt , invariant under spin-flip and conserving the total 
number of electrons and the magnetization, was consid- 
ered in some papers pd|. Its Hamiltonian reads 
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where /Ze is the chemical potential, ^, Cj^^ are fermionic 
creation and annihilation operators ({cj-.^' , Cfe,cr} = 0, 

) on a l-dimensional chain A (j, k £ A, a £ {TiJ-Dj 
and < j, k > stands for nearest neighbors (n.n.) in 
A. In (|l|) the term t represents the band energy of 
the electrons, while the subsequent terms describe their 
Coulomb interaction energy in a narrow band approxi- 
mation: U parametrizes the on-site diagonal interaction, 
V the neighboring site charge interaction, X the bond- 
charge interaction, W the exchange term, and Y the pair- 
hopping term. Moreover, additional many-body coupling 
terms have been included in agreement with pi[ : X cor- 
relates hopping with on-site occupation number, and P 
and Q describe three- and four-electron interactions. 
It is easily verified that Ti. commutes with the global sym- 
metry algebra su(2) ® it(l), where su{2) is generated by 

S+ = <^].f']A^ = '5'+, S^ = \ Ej("j,T - and 
u(l) corresponds to the conservation of the total number 
of electrons Me — X]-j("j,T + '^j.i)- 

The Hamiltonian (|l]) acts on a vector space V which is 
the tensor product of 4-dimensional on-site spaces. The 
local vector space associated with the j— th site is 

spanned by the basis vectors \oL)j — Yi-^^\o)j , a = 1, ..,4 
where /i^^^ - • /t*^' - • ~ 1 ■ h^^^ - 
and \o)j is the local vacuum (c^-^jo) = Vcr). Due 
to the fermionic nature of the creation operators, such 
vectors are endowed with an intrinsic grading e, i.e. 
e(l) = e(2) = 1 and e(3) = e(4) = 0. The basis vectors 

of are defined as |ai, a2, ... ajv) = h'f^\..h''^"^\0), 
with |0) the global vacuum {cj^\0) — Vj, cr). In order 
to work out calculations for fermionic Hamiltonians, it 
is useful to associate the basis vectors of an ordinary 
cuclidean space with the vectors \a)j of Vj, so that 
a) ^ Ca — (0, .., 1, .., 0)-^, the 1 being at the a-th po- 
sition. The Cq's inherit the grading e of the vectors \a)j. 
Moreover, we shall also introduce a matrix representation 
H for the Hamiltonian H, whose elements are defined 
through H |/3i, /32, • ■ ■ /Jw) = -ff^i'.'.'.X" a2, . . . ajv)- 

Ti is then said to be integrable if it exists a matrix R, 
acting on C* (X> C^, which i) obeys the YBE, 

R23{U - v)Ri2{u)R23{v) = Rl2{v)R23{u)Rl2{u - v) (2) 

on C'*(8'C^CE5C'', the subindices of the i?- matrix referring 
to the spaces on which R acts nontrivially; ii) has a shift 
point, meaning that i?(0) = I; and iii) is related to the 
matrix representation H through the relation 

= ^I(8)I(g)...(8) a„i?|„^(,(g)I(g)...(X)I , (3) 

implying that the first order of R in the spectral parame- 
ter u is the representation of the (two-site) Hamiltonian. 

In a recent paper Q , it was shown that when the ma- 
trix i? is a finite polynomial of degree p in u, R{u) = 



1 + + ^RSp\ in order for the YBE to be satisfied the 
coefficients in the polynomial expansion must themselves 
obey a certain number of equations. In particular the 
highest degree coefficient order R'^'p^ must satisfy 

6(P) 6(P) _ p(p) r(p) nip) (A\ 
^23 ^12 ^23 — ^12 ^23 ^12 ■ [V 

cxI , (5) 

which are realized to be the braid group relations spe- 
cialized to the symmetric case Rij = Rji. This is also 
known as the symmetric group. 

The simplest non-trivial case of a polynomial i?-matrix 
is obtained for p = 1. In this case R{u) = l + uM^\ and, 
due to (^, R'^^^ is nothing but a braid operator. The solu- 
tions of the braid group equation have been exhaustively 
investigated for a two-dimensional local space (see for in- 
stance [^). On the contrary, to our knowledge, little is 
known for the 4-dimensional case. Therefore in the next 
section we look for 16 x 16 matrix solutions of eqns. (^- 
(^ parametrized in two (a priori) different forms: i) the 
generalized permutators (to be defined in the following), 
ii) the extended Hubbard matrices. Our results show 
that the latter class of solutions is actually a subclass of 
the former, so that every integrable extended Hubbard 
model deriving from a first order polynomial i?-matrix is 
a generalized permutator. 

It is well known that, among the solutions of the sym- 
metric group relations (||)-(||), one always finds the per- 
mutation operator P, defined by P{ea Cg) Cjs) = ep ® Ca- 
In fact, this is true whatever the dimension of the repre- 
sentation is, and also holds for any graded permutation 
operator Pg, where Pg{ea e^) — {—Y'-^^'^'^'^^ep (g) Ca. 
Therefore, both P and Pg give rise to integrable models. 
In particular, the model corresponding to the i?-matrix 
having as first degree coefficient is the EKS Hamilto- 
nian m. One may wonder whether there are other solu- 
tions of (^, (|^) generalizing the structure of P, Pg. Here 
we propose the generalized permutator ones, given by op- 
erators H which either transform one product of basis 
vectors into the reversed product, or leave it unchanged. 
In the matrix form this means that there is precisely one 
non-zero entry in each column and row of H. Moreover, 
due to (||) , the non- vanishing entries of H must be equal 
to -|-1 or —1, up to an overall multiplicative constant. 
Exphcitly, Va, /? = 1, . . . 4, 

\l{ea®ep) = ea/3 sj^^ (e^ ® e/j) -|- (I-Eq/s) (e^j e^) . 

diagonal terms off-diagonal terms 

Here eafs is a discrete (0 or 1) valued function satisfying 
f-ap = i/3a and Cqq. = 1, which selects the diagonal/off- 
diagonal terms. Moreover, sj^^ — ±1 accounts for ad- 
ditional signs of the diagonal entries, while — ±1 
stands for the signs of off-diagonal terms; we shall im- 
pose = in order to make H a symmetric matrix. 

The generalized permutator H has the form given by 
matrix (|), where tr^^ = e„^s^^ and = (1 - e^^)s°^, 
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so that a'^p 



Va,/3 {i.e. actually there is only square equals +1). 



one non- vanishing element at each row/column, and its 
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Let us denote by TV the set of index values, i.e. 
N = {1,2,3,4}; the function e^/? on A/" x A/" characterizes 
the structure of any generalized permutator. Let us de- 
note by V the subset of A/" x TV of the couples (a, (i) for 
which = 1: they determine the positions 4 (a — 1) + /3 
of the non- vanishing diagonal entries; since Cq,^ = C/jq, 
whenever cr^^ 7^ we also have that cr^^ 7^ 0; in other 
words, the subset T) can always be written in the form 
V = [J^ £i X £i where the Si's are disjoint subsets of TV. 
Inserting ageneralized permutator H into the braid group 
equation it can be easily found that 11 is a solution if 
and only if a^^ = St V (a, f3) e Si x Si, where Si — ±1 
are signs. The values of the Si's can be chosen inde- 
pendently, and the remaining off-diagonal non- vanishing 
entries are also free. In the following we shall refer to 



these solutions as the symmetric group generalized per- 
mutators and denote them by lis. 

As the above definition of Ilg is independent of the di- 
mension of the representation, we conjecture that they 
are solutions of the symmetric group relations in any di- 
mension. However, their actual number depends on the 
dimension of Vj. In the present case this is four, and the 
number of different n^'s is 1440. 

Let us now consider explicitly the possibility that the 
representation of the two-site extended Hubbard Hamil- 
tonian (up to an additive constant c)is itself a solution of 
the symmetric group relations It turns out that 

actually there are 96 different possible choices of values 
of the physical parameters in (Q) satisfying (^-(||). They 
can be cast into six groups as^) shows. 





ffi(si, ...,85) 


1 H2(S1, . . . , Ss) 


-tfsCsi, S2, S3) 


-H"4(si, S2, S3) 


-H"5(si, S2, S3) 


-H"6(si, S2, S3) 


t i 


1 


1 1 


1 


1 


1 





X 


1 


1 1 


1 


1 


1 





X 


1 + S2 


1 1 + S2 


1 -1- S2 


1 + S2 


1 


1 


u 


2si 


1 2si 


4si 


4si 


2si 


-2si 


V 


Sl 


1 Sl + S4 


Sl 


Sl -1- S3 


Sl + S3 





w 


S4 





S3 











Y 


S3 


1 S3 








S2 


S2 


P 


S4 - Sl 


1 — Sl — 2S4 


S3 - 2si 


-2(si +S3) 


-(si + S3) 





Q 1 


— 2S4 + Sl + S5 


1 4S4 + Sl + S5 


4si - 2S3 


4(si -I-S3) 


Sl + S3 


Sl + S3 


M 1 


-2si 


1 -2si 


-2si 


-2si 


-2si 





c 1 


Sl 


1 Sl 


Sl 


Sl 


Sl 


Sl 



(7) 



Here Si = ±1, i = 1, . . . , 5 are arbitrary signs; the first 
two groups consist of 32 different solutions each, while 
any of the other four groups is made of 8 different cases. 
One can notice that t = X is a common feature exhibited 
by all the solutions, implying that the number of doubly 
occupied sites is a conserved quantity for those Ti. that are 
derivable from first-degree polynomial i?-matrices. This 
feature is important in that it means that in these cases 
Ti. can be diagonalized within a sector with a given num- 
ber of up and down electrons and doubly occupied sites. 



In practice, the solvability of the model in one dimension 
is not affected by having values of U and fi in Ti other 
than those reported in (Q) (see also ||l|). This is true 
also for the term proportional to the identity: indeed the 
value c — Sl exhibited by all the solutions just stems from 
requiring that the i?-matrix is a first degree polynomial. 

Remarkably, all solutions (0) do belong to the class of 
symmetric group generalized permutators Hg. In fact, it 
is easily checked that they are all those Us which pre- 
serve the spin S and the number of electrons Afe , the lat- 
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ter conditions reducing the 1440 integrable models iden- 
tified by the n^'s to precisely the 96 cases reported in 
(^. In other words, the integrable extended Hubbard 
Haniiltonians having a first degree i?-matrix are in one- 
to-one correspondence with the symmetric group gener- 
alized pcrmutators invariant under su{2) © 
One may wonder whether the 96 solutions have non- 
similar i?-matrices, i.e. whether they correspond to in- 
dependently integrable models. In the Yang-Baxter for- 
malism, two matrices R and R' are said to be similar if 
they are connected by a similarity transformation of the 
form T — A^^A, A being defined on the site space and 
(g)* being the graded product when this is the case, it 
is proven in general that, for R satisfying YBE, R' fulfills 
the YBE as well , the integrability of the latter being 
automatically deduced from the integrability of the for- 
mer. To answer the above question, we first observe that 
the traces of the various solutions can be different, as- 
suming even values from —10 to 10. Indeed, for the first 
four groups, the trace is given by the value of Q, whereas 
for the last two groups it is given by 9s3 + Si and 9s i -I- S3 
respectively. It is then obvious that solutions with dif- 
ferent trace values can not be similar, so that the above 
question needs to be answered only within a class of given 
trace value. Let us focus for instance on the models with 
vanishing trace. They are 12 and all have Q = P = 0; 8 
of them belong to the first group and are characterized by 
the Hamiltonian _ffi(si, S2, S3, si, si), whereas the other 
4 are subcases of the fourth group with H^^^si, S2, — si). 
We have verified that no pairs of their i?-matrices are 
similar, so that the 12 models with vanishing trace are 
independently integrable models. 

However, there are other kinds of unitary transforma- 
tions mapping some of the solutions ^ into other 
ones. This feature stems from the invariance proper- 
ties of the Hamiltonian which is mapped into it- 
self by means of the canonical transformations oi^a — > 
{-Yci^a, and Ci^„ -> [1 - (1 - (-)')ni^_cr]cj,o-, yielding 

n{t,x,x,u,v,w,Y,p,Q,^l) ^ ~n{t,x,x,-u,-v, 

-W, -Y, -P, -Q, -n) and H {t, t, X, U, V, W, F, P, Q, ^i) 
^ n [t, t, 2t - X, U, V, W, -y, P, Q, ^i) respectively. 

Within the various groups one can recognize some 
known integrable models, all having P = = Q, t ^ 
0. The EKS hamiltonian |l| is i7i(-l, -1, -1, -1, -1) 
which in fact corresponds to the standard graded per- 
mutator Pg, whereas i?i(l,— 1,1,1,1) is the EKS model 
with opposite grading. The AAS hamiltonian |l^] is 
nothing but iJ4(— 1, — 1, 1), whereas in the fifth group 
Hsisi, 32, Si) are particular cases of the supersymmetric 
U model whose most general form is however asso- 
ciated with a second degree polynomial i?-matrix It 
is worth stressing that even for the more standard case 
P = = Q, t ^ 0, the integrable models presented here, 
amounting to 16, are a far wider class than those which 



can be found in the literature. 

Once for a given model the i?-matrix is known, its 
spectrum is obtained within QISM by diagonalizing the 
corresponding Transfer matrix. In practice here the ex- 
act solutions, the details of which are deferred to a future 
publication, is achieved by closely following the Bethe 
Ansatz analysis for Hamiltonian acting as pcrmutators, 
first considered by Sutherland in [|l^ (see also (l) for 
its application to the EKS model). Let us just men- 
tion one physically interesting feature of the phase dia- 
gram which emerges in all the solutions characterized by 
V — ^{W ± Y), and survives also in dimension greater 
than 1, as shown in |]ll]| : it is the appearence, below a 
critical value of U, of a superconducting phase character- 
ized by 7]-pairs | |T5| , p^ . 

Finally, it is worth noticing that the solutions to equa- 
tions (^-(||) presented here are expected to be relevant 
also in the search for completely integrable models with 
higher degree polynomial _R-matrices. Indeed they pro- 
vide a class of possible coefficients for the highest 
degree term, in the form of symmetric group generalized 
permutators. Work is in progress along these lines. 
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